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Abstract 



> 

\^ . Based on a simple but effective regularization-renormalization method (RRM), 

^ . the running couphng constants (RCC) of fermions with masses in quantum electro- 

Q^ . dynamics (QED) and quantum chromodynamics (QCD) are calculated by renor- 



malization group equation (ROE). Starting at Q = {Q being the momentum 



O-i. transfer), the RCC in QED increases with the increase of Q whereas the RCCs 



for different flavors of quarks with masses in QCD are different and they increase 



^ . with the decrease of Q to reach a maximum at low Q for each flavor of quark and 

H : 

then decreases to zero at Q — > 0. The physical explanation is given. 
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I. Introduction 

In most literatures and textbooks, tlie running coupling constant (RCC) in quantum 
electrodynamics (QED) is usually given as (see, e.g. Ref. [1] and Eq. (48) below): 

where rrie is the electron mass and a = |^. Eq. (1) is very important in physics for it 
unveils the monotonically enhancing behavior of electromagnetic coupling constant a in 
accompanying the increase of momentum transfer Q between two charged particles and 
shows the existence of Landau singularity at an extremely large Q. However, in our opinion, 
there are still three aspects that can be improved in this paper, (a) Besides electron, the 
contributions of other charged leptons and quarks can not be neglected, (b) While Eq. (1) 
is scale invariant, it ignores totally the particle mass effect which is also important at low 
Q region, (c) While the normalization in Eq. (1) is inevitably made at a{Q = rrie) = a = 
(137.03599)^^, we prefer to renormalize it at the Thomson limit {Q — > 0) irrespective of 
the particle mass. 

As for quantum chromodynamics (QCD), similarly, the RCC of quark is usually expressed 
for massless quark and so is independent of the flavor of quark. For instance, it reads [2] 
(We use the Bjorken-Drell metric throughout this paper.): 

4:71 

^s{Q) = (n2i\2 T (2) 

/3o In(gVA^CD) 

where [3q = y C^ — |n/ with Ca = 3 and Uf the number of flavors of quarks. The singularity 

of Q in Eq. (2), Aqcd, is an energy scale characterizing the confinement of quarks in QCD, 

Aqcd — 200MeV experimentally. 

While Eq. (2) successfully shows the asymptotic freedom of quarks at high Q, it is not 

so satisfying at low Q region, especially for heavy quarks. The mass of c or 6 quarks, let 

alone t quark, is much higher than Aqcd- In other words, the c (or b, or t) quark does not 

exist at low Q regin beneath the threshold for creating cc (or bb, or ti) pair and the latter 



is different for different ffavor. Therefore, instead of Eq. (2), we need a new calculation of 
renormalization group equation (RGE) for RCC to discriminate different flavors of quarks. 
Evidently, it is necessary to take the mass of quark into account. 

In recent years, based on the so-called derivative renormalization method in the literature 
[3-11], proposed by Ji-feng Yang [12], a simple but effective renormalization-regularization 
method (RRM) was used by Ni et al [13-17]. It is characterized as follows. When encoun- 
tering a superficially divergent Feynman diagram integral (FDI) at one-loop level, we first 
differentiate it with respect to external momentum or mass parameter enough times until it 
becomes convergent. After performing integration with respect to internal momentum, we 
reintegrate it with respect to the parameter the same times to return to original FDI. Then 
instead of divergence, some arbitrary constants Cj {i = 1, 2, ■ ■ ■) appear in FDI, showing 
the lack of knowledge about the model at quantum field theory (QFT) level under consid- 
eration. They can only be fixed by experiments or by some other deep reasons in theory. 
Since all constants are fixed at one-loop level, all previous steps can be repeated at next loop 
expansion. The new RRM has got rid of the explicit divergence, the counterterm, the bare 
parameter and the ambiguous (arbitrary) running mass scale /i quite naturally. In section II 
we will explain this method by calculating the RCC in QED [16,18] which also serves as the 
basis of the following sections. Then in Sec. Ill the relevant formulation of RGE for RCC 
in QCD is presented. The numerical results are given at Sec. IV. The final section V will 
contain a summary and discussion. 



II. RGE calculation of RCC in QED 

As is well known, there are three kinds of Feynman diagram integral (FDI) at one-loop 
level in QED. 

1. Self-energy of electron with momentum p 
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The FDI for self-energy of electron reads (e < 0) [19-22] 



.S(p) = He)M^f^7^ , \ ^ 7- 



2 / d'k ^ .^^ 

1 (ix 



D A;2[(p — A;)^ — -m^] jo [A;^ — 2p ■ /cx -|- (p^ — 'm^)a;]^ 
A^ = 9i.vl^{i>- h + "^)7'' = -2(j^- h) + 4m. 

We first perform a shift in momentum integration: k — > K = k — xp, so that 

- i^{p) = -e^ I dx[-2{l -x) ^ + Am]! (4) 

^0 



and concentrate on the logarithmically divergent integral 

d^K 1 



(27r)4 [K^ - M2]2 
with 

A differentiation with respect to M^ is enough to get 

dl -i 1 



2 



(5) 



(6) 



(9M2 (47r)2M2' 
Thus 

carries an arbitrary contant Ci = — In/Xg- After integration with respect to the Feynman 
parameter x, one obtains 

E(j9) = A + B ^ 



a r, "H^^ im^—p^) m?' + p^ im^—p^). 
A .-^--3 5 1 + 5 In 



B = ^0n^-3-i^i-^[l + ^i-J^ln-l-5^]}. (8) 



Using the chain approximation, one can derive the modification of electron propagator as 

i i 1 iZ^ 



m 6 — ml — ^M j6 — mn 



(9) 



Z2 = {l-B)-^ ^1 + B (10) 

rriYi = ~ (m + A)(l + i?) ~ -m + Sm 

1 ~ B 

5m ~ A + wlB. (11) 

For a free electron, the mass shell condition p^ = m? leads to 

am m"^ 

dm, = — — (5 — 3 In -^). 
All' III' 

We want the parameter m in the Lagrangian still being explained as the observed mass, i.e.. 
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?72r = mobs = f^- So 5m = leads to In ^ = |, which in turn fixes the renormalization 
factor for wave function 

Z2 = l-^. (12) 



2. Photon self-energy — vacuum polarization 

n,.(g) = -H^)^Tr/A^7.^7. ^_]_^ - (13) 

Introducing the Feynman parameter x as before and performing a shift in momentum inte- 
gration: k ^ K = k — xq, we get 

n^,(g) = -4e2 f'dxih + h) (14) 

Jo 

where 

' J (27r)4 {K^-M^Y ^ ' 

with 

M"^ = m^ + q^{x'^ -x) (16) 



is quadratically divergent while 

'~y (27r)4 {K^-M^f ^ ^ 

is only logarithmically divergent like that in Eqs. (5) — (7). An elegant way for handling Ji 
is modifying M^ into 

M'^{a)=m^ + q^{x^ -x)+a (18) 

and differentiating Ji with respect to a two times. After integration with respect to K, we 
reintegrate it with respect to a two times, arriving at the limit a ^ 0: 

h = ^,{[m' + q\x'-x)]\n'^l±^^f^^ + C,} (19) 



with two arbitrary constants: Ci = — ln/x| and C2. Combining Ji and I2 together, we find 

T-r / \ 8ie^ , 2\ /^ 7 / 2 X , ■m^ + g^(x^ — x) iAe^ ^ , , 

^p^M) = TT^yli'l'' ~ 3^^<i ) / dx{x - x) In 2 llTvlS^^^^- (20) 

The continuity equation of current induced in the vacuum polarization [19] 

q^U^M) = (21) 

is ensured by the factor (g^g,y — g^vq^)- So we set C2 = 0. Consider the scattering between 
two electrons via the exchange of a photon with momentum transfer g ^ [19]. Adding the 
contribution of Il^u{q) to tree diagram amounts to modify the charge square: 

e -^ Cr = Zse 
2 2 

The choice of /is will be discussed later. The next term in expansion when g 7^ constributes 
a modification on Coulumb potential due to vacuum polarization (Uehling potential). 

3. Vertex function in QED 



aAp'.p) = i-ierj f^.i^"-jrj^„"'^jrj^-''- (23) 



For simplicity, we consider electron being on the mass shell: p^ = p' = -m^, p' — p = q, 
p ■ q = — Y- Introducing the Feynman parameter u = x + y and v = x — y, we perform a 
shift in momentum integration: 

k^'K = k — (p-\ — )u V. 



Thus 



K^ = -te'[h^^ + h] (24) 

'' = r ^" I-u ^" / ^ HK^-M^r (25) 

M2 = (m^ - ^y + ^v^ (26) 

1 /•« r d^K A 

A^, = (4 - 4m - 2M^)m^7^ + 2?(m^ - u)mq''a^^ 

u V 
-{2-2U + — - j)q\ - (2 + 2M)t;mg^ (28) 

Set K^ = K^ — M^ + M^, then 1^ = 1^ — 3^- -^3 ^^ only logarithmically divergent and can 
be treated as before to be 



'^ = Wfl * L *■ '" ' — TT^ *''' 



with /i^ an arbitrary constant. Now g^ = — Q^ < (Q^ > 0) 



— 7 7T7 'l 1 

F(cj) = m , LJ 



1-u' ./i^TT' 



Q2 
On the other hand, though there is no ultra-violet divergence in I4, it does have infrared 

divergence at -u — ;► 0. For handling it, we introduce a lower cutoff rj in the integration with 

respect to u 

+4(2 Inr^ + -)wF{w)-f^ + [-F{w) - 2]^^}. (31) 



Combining Eqs. (30) and (31) into Eq. (24), one arrives at 

A^(p',p) = -£{[ln^-- + — F(a;)]7^-(41nr7 + 5)^F(^)m27^ 

--^F{u)mq'a^,-2{2\n7^+ -)ujF{u)^,}. (32) 

When Q^ « m^, we get 



a 11 rn? a q" a 1 4 g^ 

— ( In — T + 41nr7)7u + i crnu ( — I — mr?) — r 

47r^ 2 fij " '^ 4:71 m ^ 47r^6 3 ^' m^ 



It means that the interaction of the electron with the external potential is modified 

- 67^ -^ -e[7^ + A^(p',p)]. (33) 

Besides the important term i-^^'^fiu in ^nip'^p) which emerges as the anomalous magnetic 
moment of electron, the charge modification here is expressed by a renormalization factor 

Z{^ = l + —{[2-\n — -—F{w)] + {4\nr] + 5)—^F{w) + {2\nr] + -)2wF{w)}. (34) 

The infrared term (~ In?]) is ascribed to the bremsstrahlung of soft photons [20,22] and can 
be taken care by KLN theorem [23]. We will fix /ii and rj below. 

4. Beta function at one-loop level in QED 

Adding all three FDI's at one loop level to the tree diagram, we define the renormalized 
charge as usual [2, 20-22]: 

eR = i^Zl^^e. (35) 

But the Ward-Takahashi Identity (WTI) implies that [20-22] 

Zi = Z2. (36) 

Therefore 

ftR = ^ = Zsa. (37) 



Then set p^ = irP' in Z2 and Q^ = in Zi with /ii = /i2, yielding 



lnr7 = -|. (38) 



For any value of Q, the renormalized charge reads from Eqs. (20) — (22): 

e^{Q) = e{l + - /' dx[{x - x^) In ^'^"^ ~ f ^ + "^' j} (39) 

TT JO /U3 

The observed charge is defined at Q^ ^ (Thomson scattering) limit: 

Cobs = eR|Q=o = e (41) 

which dictates that 

Ai3 = rn. (42) 

We see that e\{Q) increases with Q^. For discussing the running of cir with Q^, we define 
the Beta function: 

/9(«,g) = g^«R(Q) (43) 

From Eq. (39), one finds: 

, 2a^ 4a^m% 2m' ^Q^ j, AQ'm' - Q\ 

^^"'^^ 37r 7rg2 ^ + v/Q^T4g^^ VQ^T4g^^+g2> ^44j 

2«2 O^ O^ 

/3(a,g)^^^, (^«1) (45) 

loTT m^ 4m^ 

^. la' 4a'm' Am' . ._ . 

which leads to the well known result P{a) = ^ at one loop level at Q' -^ cxd. 

5. The RGE in QED with contributions from 9 kinds of fermions with masses 

Usually, the GRE in QED is obtained by set Q — > cxd and a — > o.r{Q) in the right 
hand side of Eq. (43), 

4"- If- '*^' 
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Then after integration, one yields analytically (see Eq. (1)): 

a 

Stt m 

However the renormalization is forced to be made ak, Q = m so that 

«i?|Q=m = a- (49) 

We are now in a position to improve the above GRE calculation in three aspects as 
indicated at the beginning of this paper. For constructing a new GRE, we replace the 
constant a in right hand side of Eq. (44) by aR{Q) and add all the contributions from 
charged leptons and quarks together, yielding: 



dQ i Stt nQ 



2m , ^Q^ + AQ^mj-Q' 
1 + , In \ 



(50) 

where 

r 1, i = e,n,T 

^i = { 3x(f)2 = |, t = u,c,t (51) 

[ 3x (-|)2 = |, i = d,s,b. 

Adding up contributions from particles with mass rrie, m^, rrir, fUc = l.OSlGeV, nib = 
4.326GeV , rrit = 175GeV we calculate the running coupling constant numerically from 
aji{Q = 0) = a till 

a^(g = mz = 9l.l884:GeV) = (131. 51)-^ (52) 

in comparision with the experimental value [24], 

aexp(Q = mz) = (128.89)-^ (53) 

The remaining discrepancy is ascribed to the contribution of light quarks {u, d, s) with 
average mass 

Thg = 92MeV, q = u,d, s. (54) 

If we adopt the following values for the mass of light quark: 

ruu = 8MeV, m^ = lOMeV, m, = 200MeV (55) 
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which are not far from the ratios found by Yan et al. [25] via the analysis of mass spectuum 
of mesons, then the fit will be rather good. See Fig. 1. 



III. RGE of RCC in QCD 

1. Self-energy of quark with mass rrii 

For convenience, we use the notation and diagram in Ref. [2] at one-loop level. Then the 
self-energy of quark with momentum p reads 

i:,{p) = -i{A, + Bi i)). (56) 

The similar procedure as in previous section leads to the renormalization constant for wave 
function: 

^2, = (l-5,)-i^l + 5,(p,m,) (57) 



«,„„„„., ra\ 



2 



{mj - p^) {mf + p^) (mf - p^) - 



Z2, = 1 + -^T'^T'^Iln^ - 3 - ^"^' " ^ [1 + ^"^' ; "M n ^'"^- / ' ]} (58) 

All IJ,2i P^ P^ ITT'i 

where a^ = |^ is the strong coupling constant, T"'T"' = |, and fi2i is an arbitrary constant 
like that in Eq. (7). 

2. Self-energy of gluon 

The combination of contributions from the gluon loop and the Faddeev-Popov ghost field 

leads to 

ids 5 Q^ 

^%,abiQ) = -r^abCA-i9f,uQ^ + Qt^qu) In -^ (59) 

^ 4:71 O fl^ 

where Q^ = —q^ > 0, Ca = 3, and /is being an another arbitrary constant (See Eq. (20)). 
The third contribution is coming from quark loop with mass rrii {i = u, d, s, c, b, t): 

^%,abiq) = ^Sabiq^q^ - g^vq^) / dx{x^ - x) In — 2 (60) 
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(the quark notation qi should not be confused with the momentum transfer q). 
Combination of Eq. (59) with (60) induces the change of as'- 

as — > Z^as 

with 

Z3 = 1 + -— --C;Aln^- + > 4 / dx{x-x)\n ^ . (61) 

3. Vertex functions in QCD 

There are two kinds of vertex function for one species of quark with mass rrii at one-loop 
level in QCD, T^^^{q) and fg^(g) (see Ref. [2]): 



jj,^ 2 2uJ' 

2UJJ _, , n „,„ , 7, 



rS'w - S(f -r-ndin^-iH-^^Mh,. 



(4 In r, + 5)-^F(L0i)mi'f„ - 2(2 In ,, + -)LO,F(u,)'t^} (62) 



fs'c) - sf r*/->"(-^^>('4+ 



^2 



+7. ^, } (63) 

where Hi {rf) is an arbitrary constant introduced for dealing with the ultraviolet (infrared) 
divergence (see Eqs. (23) — (34)), 

Ui = , ^, F{uJi) = In- (64) 

'l + 4m2/Q2 1 - ^, 



^t 



M^ = mf{l-uf + ^{u^-v^). (65) 

Here the new renormalization method has been used and two terms related to the anoma- 
lous magnetic moment of quarks have been omitted. The two Feynman diagrams give the 
correction of vertex function at one-loop level 

- igsT'^^, -^ -tgsT'^il, + r^V + r2^) = -tgsT^l,/Zu. (66) 
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Then, 



-2(21„,+ -V.F(..)} + j^f ^ *L *"{-^'"7t-2 



4. Beta function at one-loop level in QCD 

Combining all of the above one- loop Feynman diagrams and considering p = | in Z2i, 
the strong coupling constant ag is modified to 

Z'^ Z 

as — > asi{Q,mi) = ^' ag. (68) 

For discussing the running of asi{Q,fni) with Q^, we define the /3-function 

d d 

Pi{Q,mi) = Q-—^asi{Q,mi) = 2Q^-—asi{Q,mi) 

•9 „o d „ d 



By denoting 



2Q^".(aQj^l.+ aQjZ3+gQjZ,f). (69) 



^ii = -7-7^B2i{Q,mi) 



(9 „ a 



^B3{Q,mur--,'mt) (70) 



dQ^ ' 47rg2 



^li = TT^^«('5>"^i)> 



we get 



ag2 i* 47rg2 



A(g, m,) = ^{Bu + B2^ + 53). (71) 

ZTT 



5. RGE for quark qi with mass rrij in QCD 

The RGE is established by simply substituting the Ug by asi{Qimi) at the right side, 
yielding 

Q^asiiQ^nii) = l-{Bu + B2^ + B3)aUQ,m). (72) 
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IV. Numerical calculation of RGE in QCD 

Obviously, Eq. (72) can only be integrated numerically for one species of quark with 
mass rrii. We adopt the experimental data Q = rriz = 91. 18846*61^, a^j = 0.118 [26,27] as 
the initial value of integration. Then, asiiQjirii) becomes 

oils + 2;; Jq y^u + ^2i + ^3)QdQ 

where 

r^ /^ N I ,mf ^, . mf ,^ 4m? o^, , ,1 2m? n 1, 

-9 + 3 / duGi{u,Q) (74) 

Jo 



4"^f/i \/ o 2 1 M ^i + " , 1/ 2 , 4mf 



G.(«, Q) = ^(1 -«)(«- 2«2 -—) In ii^ + -^(«2 + ^«(l -«)) (75) 



^il-uy + u\ u;,= \ F(c.,)=lni±^, (76) 

\ Q^ yi + 4m2/Q2 l-o^i 



3 Q^ Q"' m^ 



52.(g,m,) = -(1 + —^(-1 + -^ln(l + ^))) (77) 



* 4m? 8m^ 
S3(Q,m.,---,m,) = -l-5:(-^-^.;,F(c..)). (7J 

i=u ^ ^ 

The results are shown in Figures 2 and 3. 



V. Summary and discussion 

1. Let us first check the zero mass limit of above equations for returning to the familiar 
result Eq. (2). For the purpose we look directly at the Zi in the limit rrii/Q — > 0, yielding 

Zri = l-f(C^ + T«T'^)ln^ 

Z2 = l + fT«TMn^ (79) 

Zl'' = ^ + mCf-ICA)lnf 
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where we have chosen In rj = —1 with another constants ;Ui = yU2 = /is = /i. This recipe 
amounts to define the value of as at high Q hmits. 
Substituting Eq. (79) into Eq. (71), we obtain 

/9(Q) = -^/3o, /3o = yCA-^n/, Ca = S. (80) 

Then the RGE reads 

Q^^RiQ) = -^PoaliQ) (81) 

with its solution precisely giving Eq. (2). 

2. Alternatively, we manage to keep the quark mass in all Bi to get the RGE (72) before 
setting the limit rrii — > 0: 

B2 — ^ 2T"T'^ 

B3 — * ~3Ca + -^nf 

Bi — > 2(^ - T'^T'') - 2Ca- 



Thus, in the limit nii — > 0, 



B, + B, + B,—. ^Uf - ^Ca = -Po- (82) 

It is interesting to compare (82) with (80), showing that 

(3o-f3o = Ca (83) 

which is stemming from the different order of taking limit: either rrii — > before the 
derivative ^2 or vice versa. 

3. But the zero mass limit is certainly not a good one as discussed in the introduction. 
And this is why one usually had to take n/ = 3 in /Sq. The mass of c or 6 quark is too heavy 
to be neglected. Therefore, we have calculated seriously the RGE for five quarks {u, d, s, c, b) 
with masses except t quark. The latter is too heavy to be created explicitly in the energy 
region considered. Notice that, however, the contribution of t quark is still existing in the 
function B^, Eq. (78). 
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4. The prominent feature of our RGE calculation is the following: 

(a) The RCC asi{Q,Tni) has a flavor dependence, i.e., it is different for different quark 
with different rrii. 

(b) The value oiasi{Q, m,i) increases from normalized value 0.118 at Q = Mz = 91.1884:GeV 
with the decrease of Q until a maximum a"^""^ is reached at Q = Aj. The smaller the rrii 
is, the smaller the Aj is and the higher the value of a™""^ will be. When Q — > 0, all a si 
approach to zero. 

(c) The value of Aj could be explaned as the existence of a critical length scale Li of qiQi 
pair 

Li ~ h/A, (84) 

while the value a"^""^ may correspond to the excitation energy for breaking the binding qiQi 
pair, i.e., the threshold energy scale against its dissociation into two bosons: 



rpthr „ max It ^ max A /*, 



(85) 



The numerical estimation of these values is listed at the table 1. It is interesting to see that 
j^thr £qj, ^^ ^ quarks is of the order of vr meson while that for c or 6 quark could be compared 
with the D^D^ or B^B^ threshold respectively. 





u 


d 


s 


c 


b 


mjc2(MeV) 


8 


10 


200 


1031 


4326 


A,(MeV) 


18.4 


18.4 


290 


1640 


7040 


max 


12.43 


9.368 


0.3027 


0.2038 


0.1610 


L.(fm) 


10.73 


10.73 


0.6809 


0.1204 


0.02805 


Ef''(MeV) 


228.7 


172.4 


87.77 


334.3 


1133 



Table 1 



Acknowledgements 

We thank Prof. Xiao-tong Song and Dr. Ji-feng Yang for discussions. This work was 
supported in part by the NSF of China. 

16 



References 

[1] I. Aitchison and A. Hey, Gauge Theories in Particle Physics (Adam Hilger LTD, Bristol, 
1982), pp 284, 287. 

[2] R.D. Field, Application of Perturbative QCD (Addison- Wesley Publishing Company, 
1989). 

[3] H. Epstein and V. Glaser, Ann. Inst. Hemi. Poincare, 19 (1973) 211. 

[4] J. Collins, Renormalization (Cambridge University Press, 1984). 

[5] G. Scharf, Finite Electrodynamics ( Springer- Ver lag, Berlin, 1989). 

[6] J. Glimm and A. Jaffe, Collective Papers, Vol.2 (1985). 

[7] M. Dtitch, F. Krahe and G. Scharf, Phys. Lett. B258 (1991), 457. 

[8] D.Z. Freedmann, K. Johnson and J.I. Lattore, Nucl. Phys. B371 (1992), 353. 

[9] P.E. Haagensen and J.I. Latorre, Phys. Lett. B 283, (1992) 293. 
[10] G. Dunne and N. Rius, ibid 293 (1992) 367. 
[11] V.A. Smirnov, Nucl. Phys. B 427 (1994) 325. 
[12] Ji-feng Yang, Thesis for PhD. (Fudan University, 1994); Preprint, (1997) |hep-th/9708T04 



Ji-feng Yang and G-j Ni, Acta Physica Sinica (Overseas Edition), 4 (1961) 88. 

[13] Guang-Jiong Ni and Su-qing Chen, Acta Physics Sinica (Overseas Edition) 7 (1998) 401; 

Internet, |hep-th/9708T55 . 



[14] G-j Ni, S-y Lou, W-f Lu and J-f Yang, Science in China (Series A) 41, (1998) 1206; 
Internet, |hep-ph/980T^U^. 



[15] G-j Ni, S-q Chen, W-f Lu, J-f Yang and S-y lou, <^Frontiers in Quantum Field Theory 3> 
Edit: C-Z Zha and K. Wu (World Scientific, 1998), pp 169-176. 

[16] G-j Ni and H. Wang, ^Physics Since Parity Symmetry Breaking^ Edit: F. Wang 
(World Scientific, 1998), pp 436-442; Preprint, Internet, |hep-th/ 9708157 . 



[17] G-j Ni, Kexue (Science) 50(3) (1998) 36-40; Internet, |quant-ph/9806009i to be published 



in a book ^Photon: Old problems in light of New ideas^ Edit: V. Dvoeglazov (Nova Pub- 



17 



lisher, 1999). 

[18] G-j Ni and H. Wang, Journal of Fudan University (Natural Science) 37(3) (1998) 304- 
305. 

[19] J.J. Sakurai, Advanced Quantum Mechanics (Addison- Wesley Publishing company, 
1967). 

[20] J.D. Bjorken and S.D. Drell, Relativistic Quantum Mechanics (McGraw-Hill Book Com- 
pany, 1964). 

[21] P. Ramond, Field Theory: A Modern Primer (The Benjamin/ Cummings Publishing 
Company, 1981). 

[22] C. Itzykson and J.B. Zuber, Quantum Field Theory (McGraw-Hill Inc., 1980). 
[23] T.D. Lee, Particle Physics and Introduction to Field Theory (Harwood Academic Pub- 
hshers, 1981). 

[24] H. Burkhardt and B. Pietrzyk, Phys. Lett. B356 (1995) 398. 
[25] D. Gao, B. Li and M. Yan, Phys. Rev. D 56 (1997) 4115. 
[26] Maria Girone and Matthias Neubert, Phys. Rev. Lett. 76 (1996) 3061. 
[27] Michael Schmelling, |hep-ex/9701002 . 



Figure Caption 

Figure 1: 

The nine curves (see from the lowest) represent respectively the contributions to the 
running electromagnetic coupling constant from 

(1) electron e only, 

(2) e and muon (yu) only, 

(3) all charged leptons e, fi and r only, 

(4) e, fi, T and c quark only. 



(5) e, /i, r, c and b quark only, 

(6) e, fi, T, c, b and t quark only, 

(7) e, /i, r, c, 6, t and -u quark only, 

(8) e, /i, r, c, b, t, u and d quark only, 

(9) all charged leptons and quarks. 

The last curve is actually coinciding with the experimental curve denoted by dot line which 
can also be fitted by assuming three light quarks {u, d, s) having average mass 92Mev/c^. 

Figure 2: 

The running strong coupling constant curves for u and d quarks. 

Figure 3: 

The running strong coupling constant curves for s, c and b quarks. 
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